Conservation & Ecology 


RAFFLES BULLETIN OF ZOOLOGY 63: 259-269 

Date of publication: 11 June 2015 

http://zoobank.Org/um:lsid:zoobank.org:pub:5E1449A4-E2BC-425F-B080-CA6C0C52DAAC 


Phylogeography of three endemic birds of Maratua Island, a potential 
archive of Bornean biogeography 

Vivien L. Chua 1 *, Quentin Phillipps 2 , Haw Chuan Lim 1 , Sabrina S. Taylor 3 , Dency F. Gawin 1 , Mustafa 
Abdul Rahman 4 , Robert G. Moyle 5 & Frederick H. Sheldon 1 

Abstract. Maratua is an oceanic island ca. 50 km off the east coast of Borneo and home to several endemic taxa 
of birds and mammals. To detennine the phylogeographic relationships of three of Maratua’s most distinctive 
avian endemics—a shama Copsychus stricklandii barbouri, a bulbul Pycnonotus atriceps hodiernus, and a monarch 
Hypothymis azurea aeria — we compared their mitochondrial ND2 sequences with those of putatively closely 
related Southeast Asian populations. We found that Maratua’s shama and bulbul are substantially differentiated 
from populations on Borneo and elsewhere in Southeast Asia, whereas Maratua’s monarch is not. The shama’s 
phylogeographic tree indicates a sister relationship between C. s. barbouri and C. s. stricklandii of northeast 
Borneo, both of which have white crowns (ND2 distance ~2%). These two taxa, in turn, are >3% divergent from 
C. malabaricus populations, which have black crowns. The relative morphological and genetic similarity of C. 
s. barbouri and C. s. stricklandii suggests that C. stricklandii had a broader distribution in eastern Borneo in the 
mid-Pleistocene and that C. s. stricklandii has more recently moved or been restricted to its current position in 
northeast Borneo. Maratua has thus acted as a biogeographic “museum”, preserving evidence of C. stricklandii s 
former distribution. The Maratuan P. atriceps is about equidistant genetically from other P. atriceps populations. 
It is distinguished by an overall grey body plumage, whereas P. atriceps elsewhere is mostly yellow (with grey 
morphs appearing only rarely). The universality of grey birds on Maratua is likely the result of a founder effect. The 
Maratuan monarch’s genetic similarity to the mainland Bornean population suggests either that it is a recent invader 
or that substantial gene flow occurs between Borneo and Maratua in this species. The genetic and morphological 
distinctiveness of the shama and bulbul are adequate to consider them both as full species. 

Keywords. Copsychus malabaricus, Copsychus stricklandii, Hypothymis azurea, phylogeography, Pycnonotus 
atriceps, rare morphs 


INTRODUCTION 

Maratua is a small oceanic island ca. 50 km east of east 
Kalimantan, Indonesian Borneo (Fig. 1). It is home to 
only 35 resident land-bird species, one monkey, and two 
rats. Among its depauperate vertebrate fauna, however, 
are some endemic taxa with substantial potential to inform 
biogeographers about the history of species distributions and 
diversification on Borneo. Maratua is close enough to Borneo 
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to attract occasional terrestrial immigrants, but isolated 
enough to act as a “museum”, preserving the morphology 
and genetic character of faunas that inhabited eastern Borneo 
long ago (Fooden, 1995). In this respect Maratua resembles, 
on a smaller scale, the Mentawai Islands off southwestern 
Sumatra (see Fig. 2B), which have preserved a surprising 
array of early Sundaic vertebrates, including an endemic 
gibbon and a proboscis-like monkey (Wilting et al., 2012). 

The first collector to visit Maratua was Harry Cusheir Raven 
in 1912-1913. He obtained specimens of all of the island’s 
terrestrial bird and mammal species for the US National 
Museum (Miller, 1913; Riley, 1930; Kellogg, 1944; Deignan, 
1959). In 1926, E. P. Mjoberg made the second collection, 
this time for Harvard University (Bangs & Peters, 1927). 
The curators who handled these early specimens were 
impressed by the distinctiveness of Maratua’s birds and 
mammals, and they named several as new species (Miller, 
1913; Bangs & Peters, 1927; Kellogg, 1944). Most of 
these names have survived to modern times as subspecific 
epithets: viz., Macacafascicularis tua (Long-tailed Macaque) 
(Fooden, 1995), Rattus tiomanicus mara (Malayan Field Rat) 
(Musser & Califia, 1982), Aerodramus salangana maratua 
(Mossy-nest Swifllet), A. fuciphagus perplexus (Edible-nest 
Swiftlet), Pycnonotus atriceps hodiernus (Black-headed 
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Fig. 1. Location of Maratua Island off the east coast of Borneo (inset). 
Bathymetric contours are, from shallowest to deepest, -40 m, -120 
m, and -200 m. Water depth data are from GEBCO One Minute 
Grid—a global one arc-minute grid; data updated in 2008 (http:// 
www.gebco.net/data_and_products/gridded_bathymetry_data/). 
Inset from www.maps-for-free.com. 

Bulbul), Copsychus stricklandii barbouri (White-crowned 
Shama), Hypothymis azurea aeria (Black-naped Monarch), 
Anthreptes malacensis mjobergi (Brown-throated Sunbird), 
and Aplonis panayensis alipodis (Asian Glossy Starling) 
(Smythies, 1999; Gill & Donsker, 2014). Early investigators 
also recognised that endemism on Maratua is due to the 
island’s permanent isolation beyond the “100 fathom” (180 
m) mark of the Sunda continental shelf (Bangs & Peters, 
1927; Riley, 1930; Kellogg, 1944). In fact, the island was 
probably never connected to mainland Borneo, not even 
during the lowest sea levels associated with glacial events 
(Voris, 2000; Bintanja et al., 2005). 

The likely long-term isolation of Maratua raises the question 
of whether its endemic populations are differentiated 
genetically, as well as morphologically, from populations 
on “mainland” Borneo and elsewhere in Southeast Asia. To 
investigate this possibility, we compared mitochondrial gene 
sequences of three species or species groups with endemic 
representatives on Maratua: Copsychus malabaricus/C. 
stricklandii, Pycnonotus atriceps, and Hypothymis azurea. 
We selected these taxa because their Maratua populations 
are morphologically distinct from those on adjacent Borneo 
and because they were the only three taxa described as full 
species in the original descriptions of the island’s birds 
(Bangs & Peters, 1927). 

Background on Maratua Island. Maratua lies ca. 50 km 
east of east Borneo in the Derawan Islands, Berau Regency, 
East Kalimantan Province, Indonesia (Fig. 1). It is a raised 
volcanic atoll of thick coralline limestone, which supports 


a relatively rich limestone-soil forest. The 2376 ha island is 
shaped like a thin horseshoe (0.3-1.4 km wide) surrounding 
a large lagoon (ca. 29.5 by 6.5 km), which is open on the 
southeastern side. The island reaches 110 m at its highest 
point (Kuenen, 1947; Tomascik et al., 1997; Becking et 
al., 2011). Kuenen (1947) believed the island developed as 
part of a barrier reef at the edge of the continental shelf. 
As the shelf gradually submerged and the coast of Borneo 
moved westward, the island became increasing isolated. At 
a later stage, the atoll was raised by tectonic activity, which 
prevented all but the southern end from being eroded away. 
Most biological research on Maratua has focused on the rich 
marine life in the surrounding coral reefs, the lagoon, and 
some anchialine pools (Ng & Tomascik, 1994; Tomascik 
et al., 1997; de Voogd et al., 2009; Becking et al., 2011). 
Since the original collections of vertebrates 100 years ago, 
little work has been done on its terrestrial fauna. 

When referring to Maratua, early writers put the island east 
of the 100 fathom (180 m) line (Bangs & Peters, 1927; 
Riley, 1930). This observation is important because 180 
m is beyond the edge of the Sunda continental shelf when 
it is exposed by eustatic sea level drops during extreme 
Pleistocene glacial events (Voris, 2000; Bintanja et al., 
2005). At this depth, Maratua is a permanently isolated 
oceanic island. However, the isolation of Maratua is not as 
extreme as implied by simple bathymetric measurement. As 
noted by Fooden (1995:90): “[Maratua] is separated from 
Borneo by a deep trench (>180 m) that is 17 km wide; the 
water gap is narrowed by the presence of a stepping-stone 
island, Kakaban (area ca. 5 km 2 ), which intervenes in the 
middle of the trench.” Indeed, as shown in Fig. 1, Kakaban 
Island almost connects Maratua to shallow water to the west. 
When sea level descended, for example, to 125 m below its 
current depth during the last glacial maximum (Bintanja et 
al., 2005), the distance from Maratua-Kakaban to Borneo was 
much reduced. This does not mean, however, that coastal 
habitats during periods of low sea level were appropriate 
for potentially dispersing forest inhabitants. 

Review of the study taxa. Maratua’s shama, bulbul, and 
monarch are each members of widespread Southeast Asian 
species or species groups (Fig. 2). Although originally 
described as full species (Bangs & Peters, 1927), all three 
were reduced from species to subspecies without comment 
by Chasen (1935) in his summary of Sundaic birds. However, 
Chasen apparently did not realise the degree of Maratua’s 
isolation, because he described this oceanic island as part 
of Borneo along with continental islands, such as the 
Natunas, the Karimatas, Labuan, the Mantananis, Banggi, 
Balambangan, Malawali, and Pulau Laut (Chasen, 1935: 
xi). Nor had he likely seen specimens of the three taxa, 
which were in American museums (Phillipps & Phillipps, 
2011: 348). 

Until recently, White-rumped Shama ( Copsychus malabaricus ) 
was a polytypic species comprising 12 subspecies extending 
from India, across Indochina to Hainan, and throughout 
Sundaland (Collar, 2005). Recently, however, two of its 
subspecies, Andaman Shama (C. m. albiventris) and White- 
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Fig. 2. Distribution (dark areas) of the three species or species 
groups whose Maratua populations were examined in this study 
(Collar, 2005; Fishpool & Tobias, 2005; Coates et al., 2006). In (B), 
the distribution of Copsychus stricklandii stricklandii in northeast 
Borneo is indicated in black. 


has black crown feathers. The white-crowned C. s. barbouri 
is larger than C. s. stricklandii, and its tail is entirely black 
(Bangs & Peters, 1927; Collar, 2004). In phylogeographic 
comparisons of C. malabaricus/stricklandii populations in 
this paper, we used the closely related C. niger of Palawan 
as outgroup (Lim et al., 2010). 

Black-headed Bulbul ( Pycnonotus atriceps ) ranges from 
eastern India through Indochina and across all of Sundaland 
to Palawan and Bali. It is divided into four subspecies. P. a. 
atriceps occurs on the Asian mainland and large continental 
islands, whereas the other subspecies are restricted to smaller 
islands: P. a. hyperemnus on Simeulue, Nias, and Mentawai 
of western Sumatra, P. a. baweanus on Bawean Island in the 
Java Sea, and P. a. hodiernus on Maratua (Fishpool & Tobias, 
2005). As in C. malabaricus, the Andaman population was 
previously considered a subspecies (Dickinson, 2003), but is 
now recognised as a full species, P. fuscoflavescens (Fishpool 
& Tobias, 2005; Gill & Donsker, 2014). The Maratuan 
subspecies has the shiny black head and distinctive blue iris 
of P. a. atriceps, but instead of a yellow body, its plumage 
is entirely a mixture of dark grey and black. Maratua birds 
are also larger than P. a. atriceps (Bangs & Peters, 1927). 
Elsewhere in P. atriceps range, grey plumaged individuals 
occur as rare morphs (Dickinson & Dekker, 2002; Wells, 
2007; pers. comm.), except on Bawean Island in the Java 
Sea, where grey individuals are more common than yellow 
(Hoogerwerf, 1966; van Balen, 1986; Dickinson & Dekker, 
2002). The sister group of P. atriceps is a clade consisting 
of the Puff-backed Bulbul (P. eutilotus ) and Black-and-white 
Bulbul (P. melanoleucos) (Moyle & Marks, 2006). We use 
P. eutilotus as outgroup for phylogeographic comparisons 
in this paper. 

Black-naped Monarch ( Hypothymis azurea) is distributed 
from India to Taiwan and south to the Lesser Sundas. It is 
a highly polymorphic species with 23 subspecies (Gill & 
Donsker, 2014) that vary mainly in shade and extent of blue 
feathering and degree of crest and breast markings (Coates 
et al., 2006). The Maratua population (H. a. aeria ) differs 
from all other subspecies by lacking a distinctive black 
occipital patch and by possessing a narrower, less velvety, 
black crescent on its throat (Bangs & Peters, 1927; Rand, 
1970). The closest relatives of H. azurea are Sulawesi’s H. 
puella (Sulawesi) and the Philippine’s H. coelestis and H. 
helenae (Fabre et al., 2012). We use the latter two species 
as outgroups in this study. 

MATERIAL AND METHODS 


crowned Shama (C. m. stricklandii), were raised to full 
species, with C. stricklandii comprising subspecies C. s. 
stricklandii of northeastern Borneo and C. s. barbouri of 
Maratua (Gill & Donsker, 2014). On Borneo, C. malabaricus 
occupies most of the island, but is replaced in the northeast 
(essentially Sabah) and on Maratua by Copsychus stricklandii 
(Fig. 3). Because of their similarity and parapatric distribution 
(Collar, 2004; Gawin, 2014), we refer to C. malabaricus 
and C. stricklandii as a “species group” in this paper. The 
two taxa differ most noticeably in crown feathering; C. 
stricklandii has white crown feathers, whereas C. malabaricus 


We obtained DNA of Copsychus stricklandii barbouri, 
Pycnonotus atriceps hodiernus, Hypothymis azurea aeria 
from toe pads of Mjoberg’s study skins: two individuals of 
each from the Museum of Comparative Zoology, Harvard 
University (Table 1). We obtained DNA of other taxa from 
preserved tissues provided by tissues collections or by 
downloading DNA sequences from GenBank (Table 1). In 
all (Table 1), our comparisons included: for shamas—16 
individuals of Copsychus malabaricus (three subspecies), 
eight C. stricklandii (two subspecies), and one C. niger 
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Table 1. Specimens compared in this study. 


Species 1 

Subspecies 

Downloaded GenBank or Specimen Numbers 2 

Locality 

Copsychus malabaricus 

interpositus 

HQ010913.1 

Vietnam: Quang Nam 

Copsychus malabaricus 

interpositus 

HQ010914.1 

Myanmar: Sagaing 

Copsychus malabaricus 

interpositus 

HQ010921.1 

Myanmar: Sagaing 

Copsychus malabaricus 

interpositus 

HQO 10923.1 

Myanmar: Sagaing 

Copsychus malabaricus 

interpositus 

HQO 10925.1 

Myanmar: Sagaing 

Copsychus malabaricus 

suavis 

HQ010908.1 

Malaysia: Sarawak 

Copsychus malabaricus 

suavis 

HM120190.1 

Malaysia: Sarawak 

Copsychus malabaricus 

suavis 

HQO 10909.1 

Malaysia: Sarawak 

Copsychus malabaricus 

suavis 

LSUMNS B58587 

Malaysia: Sarawak 

Copsychus malabaricus 

suavis 

LSUMNS B73734 

Malaysia: Sarawak 

Copsychus malabaricus 

suavis 

HQ010911.1 

Malaysia: Sarawak 

Copsychus malabaricus 

tricolor 

HQ010912.1 

Malaysia: Terengganu 

Copsychus malabaricus 

tricolor 

HQ010933.1 

Malaysia: Terengganu 

Copsychus malabaricus 

tricolor 

HM120189.1 

Malaysia: Johor 

Copsychus malabaricus 

tricolor 

HQO 10934.1 

Malaysia: Johor 

Copsychus malabaricus 

tricolor (?) 

HQ010910.1 

Singapore 

Copsychus stricklandii 

stricklandii 

HQ010915.1 

Malaysia: Sabah 

Copsychus stricklandii 

stricklandii 

HQ010917.1 

Malaysia: Sabah 

Copsychus stricklandii 

stricklandii 

HQO 10920.1 

Malaysia: Sabah 

Copsychus stricklandii 

stricklandii 

LSUMNS B73450 

Malaysia: Sabah 

Copsychus stricklandii 

stricklandii 

LSUMNS B73454 

Malaysia: Sabah 

Copsychus stricklandii 

stricklandii 

LSUMNS B73488 

Malaysia: Sabah 

Copsychus stricklandii 

barbouri 

MCZ 235961 

Indonesia: Maratua Is. 

Copsychus stricklandii 

barbouri 

MCZ 235968 

Indonesia: Maratua Is. 

Copsychus niger* 


HM120197.1 

Philippines: Palawan 

Hypothymis azurea 

aeria 

MCZ 235920 

Indonesia: Maratua Is. 

Hypothymis azurea 

aeria 

MCZ 235927 

Indonesia: Maratua Is. 

Hypothymis azurea 

azurea 

HQ010959.1 

Philippines: Luzon 

Hypothymis azurea 

azurea 

KUNHM 13983 

Philippines: Camiguin Sur 

Hypothymis azurea 

azurea 

HQ010958.1 

Philippines: Mindanao 

Hypothymis azurea 

azurea 

HQ010956.1 

Philippines: Palawan 

Hypothymis azurea 

azurea 

HQ010957.1 

Philippines: Palawan 

Hypothymis azurea 

prophata 

HQO 1942.1 

Malaysia: Terengganu 

Hypothymis azurea 

prophata 

HQ010961.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQO 10943.1 

Malaysia: Sabah 

Hypothymis azurea 

prophata 

HQO 10944.1 

Malaysia: Sabah 

Hypothymis azurea 

prophata 

HQO 10945.1 

Malaysia: Sabah 

Hypothymis azurea 

prophata 

HQO 10946.1 

Malaysia: Sabah 

Hypothymis azurea 

prophata 

HQO 10947.1 

Malaysia: Sabah 

Hypothymis azurea 

prophata 

HQ010939.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQO 10940.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQ010941.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQ010950.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQ010951.1 

Malaysia: Sarawak 

Hypothymis azurea 

prophata 

HQO 10948.1 

Malaysia: Johor 

Hypothymis azurea 

prophata 

HQO 10960.1 

Indonesia: Sumatra 

Hypothymis azurea 

styani 

HQ010938.1 

Vietnam: Quang Nam 

Hypothymis azurea 

styani 

KUNHM 10293 

China: Guangxi 

Hypothymis azurea 

styani 

HQO 10949.1 

Myanmar: Sagaing 

Hypothymis azurea 

styani 

HQ010952.1 

Myanmar: Sagaing 

Hypothymis azurea 

styani 

HQ010953.1 

Myanmar: Sagaing 

Hypothymis azurea 3 


EF052695 

Thailand 

Hypothymis coelestis* 

coelestis 

KUNHM 19062 

Philippines: Mindanao 

Hypothymis coelestis* 

coelestis 

KUNHM 19063 

Philippines: Mindanao 

Hypothymis helenae* 

personata 

KUNHM 12536 

Philippines: Camiguin 
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Species 1 

Subspecies 

Downloaded GenBank or Specimen Numbers 2 

Locality 

Hypothymis helenae* 

helenae 

KUNHM 20342 

Philippines: Luzon 

Pycnonotus atriceps 

atriceps 

LSUMNS B36320 

Malaysia: Sabah 

Pycnonotus atriceps 

atriceps 

LSUMNS B57464 

Malaysia: Sabah 

Pycnonotus atriceps 

atriceps 

LSUMNS B58500 

Malaysia: Sarawak 

Pycnonotus atriceps 

atriceps 

KUNHM 12317 

Malaysia: Sarawak 

Pycnonotus atriceps 

atriceps 

KUNHM 12633 

Philippines: Palawan 

Pycnonotus atriceps 

atriceps 

KUNHM 12639 

Philippines: Palawan 

Pycnonotus atriceps 

atriceps 

KUNHM 12642 

Philippines: Palawan 

Pycnonotus atriceps 

atriceps 

GU112683.1 

Philippines: Palawan 

Pycnonotus atriceps 

hodiemus 

MCZ 235903 

Indonesia: Maratua Is. 

Pycnonotus atriceps 

hodiernus 

MCZ 235907 

Indonesia: Maratua Is. 

Pycnonotus eutilotus * 

monotypic 

DQ402236.1 

Malaysia: Sabah 


'Names follow the classification of Gill and Donsker (2014). *Outgroup taxa. 

2 “Downloaded GenBank numbers” refers to sequences from previous studies: DQ (Moyle & Marks, 2006), EF (Fuchs et al., 2007), GU (Oliveros & Moyle, 
2010), HM (Lim et ah, 2010), and HQ (Lim et ah, 2011). “Specimen numbers” refers to material (preserved tissues or museum skins) used to produce 
new sequences in this study. GenBank numbers of sequences produced by this study are KP208752-KP208775. Museum abbreviations: AMNH, American 
Museum of Natural History, New York, New York; KUNHM University of Kansas Natural History Museum, Lawrence, Kansas; LSUMNS, Louisiana State 
University Museum of Natural Science, Baton Rouge, Louisiana; MCZ, Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts; 
NMNH, US National Museum of Natural History, Washington, D.C.; UWBM, University of Washington Burke Museum, Seattle, Washington. 

3 Exact collecting locality and subspecies not specified (Fuchs et ah, 2007). 
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Fig. 3. ND2 phylogeographic tree and haplotype network of 
populations of White-rumped Shama, Copsychus malabaricus 
(top figure), and White-crowned Shama, C. stricklandii (bottom 
two figures). A, Maximum likelihood bootstrap tree. Numbers 
by branches are bootstrap support (70% and above). Labels are 
specimen numbers; B, Flaplotype network. Circle size indicates 
proportion of individuals with a given haplotype, and cross hatches 
indicate nucleotide changes (larger numbers of changes are noted 
in square boxes). 


Fig 4. ND2 phylogeographic tree and haplotype network of Black¬ 
headed Bulbul (Pycnonotus atriceps ) populations. A, Maximum 
likelihood bootstrap tree. Numbers by branches are bootstrap support 
(70% and above). Labels are specimen numbers. B, Flaplotype 
network. Circle size indicates proportion of individuals with a given 
haplotype, and cross hatches indicate nucleotide changes (larger 
numbers of changes are noted in square boxes). 
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Hypothymis azurea 


I Mainland 
Sabah 
Philippines 


Fig. 5. ND2 phylogeographic tree and haplotype network of Black- 
naped Monarch ( Hypothymis azurea) populations. A, Maximum 
likelihood bootstrap tree. Numbers by branches are bootstrap support 
(70% and above). B, Flaplotype network. Circle size indicates 
proportion of individuals with a given haplotype, and cross hatches 
indicate nucleotide changes (larger numbers of changes are noted 
in square boxes). 


as outgroup; for bulbuls—10 Pycnonotus atriceps (two 
subspecies) and one P. eutilotus as outgroup; and for 
monarchs—27 Hypothymis azurea (at least four subspecies) 
and two each of H. coelestis and H. helenae as outgroups. 

DNA was extracted from preserved tissues using a Qiagen 
DNeasy kit following the manufacturer’s protocol. To 
extract toe-pad DNA, we used a Qiagen QIAmp DNA 
kit and applied the methods of Moyle et al. (2011). We 
sequenced mtDNA NADH dehydrogenase subunit 2 (ND2) 
using alternate methods depending on the DNA source. 
For amplification of DNA from preserved tissues, we used 
primers L5215 (Flackett, 1996) and FITrpC (STRI), as per 
Lim et al. (2010). DNA was amplified in 25 ftL reactions 
using Qiagen HotStarTaq plus DNA Polymerase kit. For 
amplification toe-pad DNA, we designed five pairs of internal 
primers from the preserved-tissue sequences (Table 2). These 
primers were used to amplify 250-300 nucleotide sections 
of the ND2 gene, which were then concatenated. All post- 
PCR products were sent to Beckman Coulter Genomics 
for sequencing. Complementary strands were aligned using 
Geneious version R6 by Biomatters (www.geneious.com). 
GenBank numbers of tissues sequenced in this study are 
KP208752-KP208775. 

Population genetic statistics were computed using Arlequin 
3.5 (Excoffier & Lischer, 2010). Phylogeographic trees of 


Table 2. Primers used to amplify and sequence DNA from study specimen toe pads. 


Name 

Species 

Primer Sequence 

R1 

Pycnonotus atriceps 

GCAATGGCTGAGGTCAGGAT 

F2 

Pycnonotus atriceps 

TAACCCACCCAGTGTCATGC 

R2 

Pycnonotus atriceps 

CCTAGGGCGGTGGATAGGAT 

F3 

Pycnonotus atriceps 

CCCATCACTACTTACTTGCATGGC 

R3 

Pycnonotus atriceps 

GGTAGGAATCCAGTTAATGGGG 

F4 

Pycnonotus atriceps 

CTCCTACTAACCCTACTCTCCC 

R4 

Pycnonotus atriceps 

GGTTAGAGGTGTGTGGTGGG 

F5 

Pycnonotus atriceps 

CGCCTTGCATACTGCGCAAC 

R1 

Copsychus malabaricus/stricklandii 

AGTGGAATGGGGCTAGTCCT 

F2 

Copsychus malabaricus/stricklandii 

AACCTGTCCTACAGCATGCC 

R2 

Copsychus malabaricus/stricklandii 

ATCCTCCTAGGGCTGTGGAG 

F3 

Copsychus malabaricus/stricklandii 

CCCCACACTACTAGTTACCATAGC 

R3 

Copsychus malabaricus/stricklandii 

TTAAGGGCGAGGAATACGGC 

F4 

Copsychus malabaricus/stricklandii 

CCTACAGCCCCAAATTAGCC 

R4 

Copsychus malabaricus/stricklandii 

ATTGTTGCTGCTGGGGCTAT 

F5 

Copsychus malabaricus/stricklandii 

AC AGGCTT CCT CCCC AAAT G 

R1 

Hypothymis azurea 

TGAGCCTTGCAGGACTTCTG 

F2 

Hypothymis azurea 

TGCCTGATCCTAACAGCAGC 

R2 

Hypothymis azurea 

CTATTCATCCTCCTAGGGCAG 

F3 

Hypothymis azurea 

CCTCGCACTCACTCAACCCAAC 

R3 

Hypothymis azurea 

GAGAAATACGGCTGCAGTGATTAGGG 

F4 

Hypothymis azurea 

CTACAACCCAAAGCTCACCC 

R4 

Hypothymis azurea 

GCCTAGCAGGGATAGCAGTG 

F5 

Hypothymis azurea 

CAGGACATAGCCCCTACAGC 
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each species or species group were estimated by maximum 
likelihood (ML) using Garli 2.0 (Zwickl, 2006). Nodal 
support was tested with 200 bootstrap replicates. Appropriate 
ML substitution models were selected based on the Akaike 
Information Criterion implemented in jModelTest 2 
(Guindon & Gascuel, 2003; Darriba et ah, 2012). The best- 
fit substitution model for P. atriceps was TIM2+G, for C. 
stricklandii, GTR+G, and for H. azurea, TrN+I. Bootstrap 
consensus trees and bootstrap values were calculated using 
SUMTREES (Sukumaran & Holder, 2010). Uncorrected 
proportional genetic distances among populations within 
each species or species group were generated using Mega5 
using both complete and pairwise divergence (Tamura et 
ah, 2011; Fregin et ah, 2012). Median joining haplotype 
networks were constructed with the program Network v. 4.6 
(Bandelt et ah, 1999) to visualise geographic relationships 
within each species. 

RESULTS 

ND2 sequence lengths of individuals in the three species 
ranged from 841 to 1,041 nucleotides. Copsychus 
malabaricuslstricklandii had 65 (6.2%) variable sites and 51 
(4.9%) parsimony-informative sites, Pycnonotus atriceps 192 
(18.4%) variable sites and 25 (2.4%) parsimony-informative 
sites, and Hypothymis azurea 64 (6.1%) variable sites and 29 
(2.8%) parsimony-informative sites. We found 18 haplotypes 
in 24 individuals of C. malabaricus/stricklandii, 9 haplotypes 
in 10 individuals of P. atriceps, and 26 haplotypes in 27 
individuals of H. azurea (Table 3). Average divergence 
values (Table 4) between Maratua and the next closest 
populations were: Copsychus stricklandii barbouri to C. s. 
stricklandii (Sabah) 2%, Pycnonotus atriceps hodiernus to 
P. a. atriceps (Sabah) 2.5%, Hypothymis azurea aeria to H. 
a. prophata (Sabah) 0.6%. 

ML trees and haplotype networks depicted three different 
phylogeographic patterns in the three species. The Copsychus 
malabaricus/stricklandii tree (Fig. 3A) had the most 
hierarchical structure, with the Maratua population (C. s. 
barbouri) as sister to the Sabah population (C. s. stricklandii), 
and this clade as distinct from C. malabaricus. Within 
C. malabaricus, the Sarawak population (C. m. suavis) 
was distinct from the Malay Peninsula and Indochinese 
populations. The Copsychus malabaricus/stricklandii 
haplotype network largely mirrored the phylogeographic 
tree (Fig. 3B). One individual from the mainland (red color 
in Fig. 3B) is nine steps removed from the other mainland 
taxa. This individual (Table 1; GenBank no. HQ010910.1) 
is from Singapore (no specimen) and possibly had been 
imported from elsewhere in SE Asia. The Pycnonotus 
atriceps tree and network (Fig. 5) indicated that Maratua 
individuals are well differentiated from members of other 
populations, whereas the other populations showed little 
phylogeographic structure. Three of the Philippine individuals 
formed a marginally distinct clade, while a fourth Philippine 
individual was largely differentiated from the others. The 
Maratua population of Hypothymis azurea was not distinct 
phylogeographically from the other sampled populations 
(Fig. 6). Although several branches leading to small groups 



Fig. 6. Borneo map showing possible Pleistocene rainforest refuges 
(yellow and orange) and the modern distribution of Copsychus 
stricklandii (orange). 


of individuals had moderate support, there was no clear 
hierarchical structure in the tree or network. 

DISCUSSION 

Biogeography and plumage patterns. All three of our 

study taxa, Copsychus malabaricus/stricklandii, Pycnonotus 
atriceps, and Hypothymis azurea, are widely distributed in 
Southeast Asia and have morphologically distinct populations 
on isolated islands. It is tempting, therefore, to suggest that 
some of the plumages found on oceanic islands such as 
Maratua may represent morphology of an early wave of 
Sundaic colonisers, whereas those on continental islands 
represent later arriving colonisers. This possibility was 
raised by Fooden (1995), who suggested that dark pelaged 
Long-tailed Macaques ( Macaco fascicularis) that occur on 
oceanic islands surrounding Sundaland, including Maratua 
(M. f. tua), are relicts of early inhabitants of the region. 
These dark forms were presumed to have been replaced on 
all continental islands by M.f fascicularis, a lighter pelaged, 
possibly more recent invader (also see Abegg & Thierry, 
2002). However, plumage variation across our three study 
taxa is complex, and in only one case, the shama, does the 
plumage pattern occurring on Maratua provide insight into 
a possible ancestral distribution. Morphology and genetic 
variation in the Maratuan monarch and bulbul, although 
perhaps not indicating ancestral distributions, still provide 
logical explanations for the evolution of their plumages on 
Maratua. 

We suggest the following scenario for the biogeographic 
history of Copsychus stricklandii populations. The white 
crown of C. stricklandii appears to have developed in an 
eastern Bornean rainforest refuge in the early Pleistocene. 
At the time of this development, C. stricklandii' s ancestors 
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Table 3. Sampling and diversity values for populations: n = number of individuals, h = number of haplotypes, var = number of variable 
sites, Hd = haplotype diversity, SD = standard deviation and Jt = nucleotide diversity. 


Population 

n 

h 

var 

Hd ± SD 

Jt ± SD 

Pycnonotus atriceps 

Maratua 

2 

1 

0 

0 

0 

Sarawak 

2 

2 

4 

1.0 ± 0.50 

0.0038 ± 0.0043 

Sabah 

2 

2 

4 

1.0 ± 0.50 

0.0039 ± 0.0043 

Philippines 

4 

4 

18 

1.0 ± 0.18 

0.0086 ± 0.0060 

mean 

2.5 

2.25 

6.5 

0.75 

0.004075 

Hypothymis azurea 

Maratua 

2 

1 

0 

0 

0 

Sabah 

5 

5 

17 

1.0 ± 0.13 

0.0078 ± 0.0051 

Sarawak 

6 

6 

20 

1.0 ± 0.10 

0.0072 ± 0.0045 

Mainland 

9 

9 

25 

1.0 ± 0.05 

0.0314 ± 0.0172 

Philippines 

5 

5 

18 

1.0 ± 0.13 

0.0340 ± 0.0210 

mean 

5.4 

5.2 

16 

0.8 

0.01608 

Copsychus malabaricus / stricklandii 

Maratua 2 

1 

0 

0 

0.0050 ± 0.0033 

Sabah 

6 

6 

9 

1.0 ± 0.10 

0.0012 ± 0.0011 

Sarawak 

6 

3 

3 

0.93 ± 0.12 

0.0057 ± 0.0034 

Mainland 

10 

8 

11 

0.96 ± 0.06 

0.0050 ± 0.0033 

mean 

6 

4.5 

5.75 

0.72 

0.003 


would have been allopatric relative to the black-crowned C. 
malabaricus, which was most likely restricted to western 
Sundaland (and perhaps western Borneo) in the early to mid- 
Pleistocene. This scenario is based on a variety of evidence. 
The genetic distance between the white- and black-crowned 
shamas is substantial (>3%; Table 3; Lim et al., 2010; Lim 
et al., 2011), especially in view of the overall similarity 
of these taxa (Collar, 2004). Also, the two shamas exhibit 
limited hybridisation where they currently come in contact 
near the Sabah border (Gawin, 2014). (This is known for the 
Malaysian side of Sabah; hybridisation has not been studied 
on the Indonesian side.) Furthermore, allopatric divergence 
in the two shama populations accords well, spatially and 
temporally, with phylogeographic evidence of Pleistocene 
subdivisions in several other groups of vertebrates caused 
by presumed dry habitat in central Sundaland (Bird et al., 
2005; Wilting et al., 2012). These other groups include not 
only birds (Lim et al., 2011), but also reptiles, amphibians, 
and several orders of mammals (e.g., Brandon-Jones, 1998; 
Gorog et al., 2004; Meijaard & Groves, 2004; Wilting et 
al., 2012). We also know that the rainforest required by C. 
stricklandii existed continuously in the mountains and some 
coastal areas of eastern Borneo since the Miocene, even 
when other parts of Sundaland were drier (Morley, 2012; 
Morley & Morley, 2013). The locations of refuges for C. 
malabaricus in western Sundaland, however, are less well 
established (Fig. 6). The Mentawai Islands and mountains 
of Sumatra are the most likely places for refugial rainforest 
because of orographic precipitation in those locations (e.g., 
Brandon-Jones, 1998; Wilting et al., 2012). Among possible 
western refuges on Borneo, the Schwaner Mountains and 
associated lowlands were suggested for lowland rodents 
(Gorog et al., 2004). 


In respect to Maratua, the mid-Pleistocene refuge of white- 
crowned shamas was likely near the Bornean coast adjacent 
to the island, and it would have provided the stock for the 
Maratua’s invasion (Fig. 6). The separation of C. s. barbouri 
(Maratua) from C. s. stricklandii (mainland Borneo) occurred 
early enough for the two populations to accrue a genetic 
distance of ca. 2% (Table 3). Currently, C. s. stricklandii 
is restricted to northeastern Borneo (essentially Sabah). Its 
positional shift was probably in response to a redistribution 
of refugial rainforest, as evidenced by the co-occurrence of 
many vertebrate rainforest endemics in Sabah (Sheldon et 
al., 2001; Sheldon et al., 2009; Gawin et al., 2014). More 
recently, presumably following expansion of rainforest in 
central Sundaland (Cannon et al., 2009), C. malabaricus 
moved from western Borneo/Sundaland across Borneo to 
meet C. stricklandii , leading to their parapatric distribution 
roughly in alignment with Sabah’s border. 

Pynonotus atriceps hodiernus, the other genetically 
distinct Maratuan taxon, appears to have had a different 
biogeographic history than the shama, at least in respect to 
its plumage evolution. Individuals with grey plumage like 
those on Maratua appear as rare morphs in various parts 
of the species’ range, although none has been recorded 
on Borneo (Smythies, 1999; Dickinson & Dekker, 2002; 
Wells, 2007; pers. comm.). Apart from Maratua, the only 
location in Sundaland where grey individuals are common is 
Bawean Island in the Java Sea (Fig. 3B; Hoogerwerf, 1966; 
van Balen, 1986; Dickinson & Dekker, 2002). Why dark 
morphs exist in Sundaland has always puzzled ornithologists. 
They occur not only in P. atriceps, but also in some other 
Sundaic rainforest passerines (e.g., Terpsiphone paradisi 
and Philentoma pyrhoptera). Studies of yellow versus black 
plumages in bananquits ( Coereba flaveola), which appear 
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Table 4. Uncorrected proportional ND2 divergence values (complete deletion of missing nucleotides/pairwise deletion of missing nucleotides). 


Locality 

Maratua 

Sabah 

Sarawak 

Palawan 

Copsychus malabaricus/stricklandii 
Sabah 

Sarawak 

Malay Pen. 

0.02/0.02 

0.036/0.036 

0.032/0.032 

0.029/0.030 

0.023/0.023 

0.02/0.02 


Pycnonotus atriceps 

Sabah 

Sarawak 

Palawan 

0.023/0.023 

0.025/0.025 

0.026/0.026 

0.005/0.005 

0.008/0.008 

0.009/0.009 


Hypothymis azurea 

Sabah 

Sarawak 

Palawan 

Mainland 

0.007/0.008 

0.008/0.008 

0.009/0.008 

0.010/0.009 

0.007/0.008 

0.011/0.011 

0.012/0.012 

0.012/0.011 

0.012/0.011 

0.013/0.012 


on different islands in the Caribbean, have shown that 
the difference is based on a single nucleotide substitution 
influencing melanin synthesis (Theron et al., 2001). In other 
cases, grey and yellow plumages are directly related by 
relative levels of lutein concentration (Brush & Johnson, 
1976; Remsen & Graves, 1995). Thus, plumage change 
in P. atriceps is probably the result of a relatively simple 
recessive genetic difference, which is expressed on Maratua 
and Bawean because of founder effect and drift. In larger 
P. atriceps populations, the recessive grey morph would be 
swamped by the dominant yellow plumage. 

We have not examined the phylogeography of the Bawean 
population of P. atriceps, but Bawean has a different history 
than Maratua in respect to isolation and habitat. Unlike 
Maratua, Bawean is a continental island that would have 
been attached to Borneo and Java during multiple glacial 
events, most recently during the last glacial maximum (ca. 
18 ka). As such, its bulbul population would have been 
subject to substantial gene flow. However, unlike eastern 
Borneo, Bawean lies in a region that would have been drier 
and probably devoid of rainforest at various times in the 
Pleistocene (Bird et al., 2005; Morley, 2012). Thus, its bulbul 
population may have been subject to multiple bottlenecks, 
invasions, drift, and other small population effects, allowing 
gray plumage to predominate. 

Maratua’s Hypothymis azurea population, although 
possessing a somewhat different plumage than the population 
on adjacent Borneo, is genetically similar not only to its 
Bornean neighbors but to all other populations we examined 
(Table 4, Fig. 6). Its lack of genetic distinctiveness would 
seem to reflect//, azurea’ s substantial dispersal capabilities. 
Its morphological distinctiveness on Maratua and elsewhere, 
resulting in 23 morphologically defined subspecies, would 
suggest a highly plastic plumage. 

Taxonomy. We have not evaluated the degree of gene flow 
influencing Maratua’s populations of shama and bulbul, which 


would provide a better assessment of their status as biological 
species than our mitochondrial dataset (Lim et al., 2014). 
Nevertheless, the genetic and morphological distinctiveness 
of Copsychus stricklandii barbouri and Pycnonotus atriceps 
hodiernus on Maratua are enough to call them phylogenetic 
species, and probably even biological species. 

Conservation. The Maratuan Shama has become noticeably 
rarer in the last few years. The Indonesian government 
is developing the island as a tourist destination and has 
recently built an airport. The influx of construction workers 
has apparently resulted in the capture of many birds for the 
pet trade, causing what was once an abundant forest species 
to become rare and in danger of extirpation (QP, personal 
observation). Even without this immediate danger, Maratua’s 
shama and bulbul deserve special protection because of their 
genetic distinctiveness, inherently small population sizes, 
and interesting evolutionary trajectories. 
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